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ABSTRACT 

We study the cosmological information contained in the Minkowski Functionals 
(MFs) of weak gravitational lensing convergence maps. We show that the MFs pro- 
vide strong constraints on the local type primordial non-Gaussianity parameter /nl- 
We run a set of cosmological A^-body simulations and perform ray-tracing simulations 
of weak lensing, to generate 100 independent convergence maps of 25 deg^ field-of-view 
for /nl = —100, and 100. We perform a Fisher analysis to study the degeneracy among 
other cosmological parameters such as the dark energy equation of state parameter w 
and the fluctuation amplitude erg. We use fully nonlinear covariance matrices evaluated 
from 1000 ray-tracing simulations. For the upcoming wide-field observations such as 
Subaru Hyper Suprime-Cam survey with the proposed survey area of 1500 deg2, the 
primordial non-Gaussianity can be constrained with a level of /nl ~ 80 and w ~ 0.036 
by weak lensing MFs. If simply scaled by the effective survey area, a 20000 deg2 lensing 
survey using Large Synoptic Survey Telescope will give constraints of /nl ~ 25 and 
w ~ 0.013. We show that these constraints can be further improved by a tomographic 
method using source galaxies in multiple redshift bins. 



Subject headings: Cosmology: cosmological parameters — large-scale structure of Uni- 
verse 
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INTRODUCTION 



An array of rece nt precise cosmologica l observations such as the cosmic microwave background 



(CMB) anisotrop ies ( Komatsu et alj[201ll ) and the large-scale structure (e.g. Tegmark et al. 



2006 



Reid et al.ll2010bl ) established the so-called standard cosmological model, in which the energy con- 
tent of the present-day universe is dominated by dark energy and dark matter. The standard 
model assumes that the primordial density fluctuations are generated via inflation in the very early 
universe, which seeded eventually all the rich structures of the universe we see today. The den- 
sity fluctuations generated through the mechanism generally follow Gaussian statistics. Deviations 
from primordial Gaussianity would provide interesting information on the early universe. For ex- 
ample, some inflationary models predict gener ation of non-Gaussi an density fluctuations, called the 
primordial non- Gaussianity (for a review, see lBartolo et al.ll2004l ). 



For the so-called local type non-Gaussian models, the initial curvature fluctuation $ is ex- 
pressed by the Taylor expansion of a Gaussian field <I>g 



$(x) = ^g{x) + /nl {^l{x) - {^l{x))) + ■■■. 



(1) 



The leading coefficient /nl determines the strength of the non-Gaussianity (iKomatsu &: Spergel 
200ll ). The amplitude of /nl depends on the perturbation generation mechanism and hence on the 



physics of inflation. 

There have been several observational probes of /nt,. For example, the b ispectrum of the 



CMB anisotropies is shown to be a powerful probe of /nl (jKomatsu et al.ll201ll ). The abundance 
and the distribution of galaxies and galaxy clusters can be also used to constrain /nl- Primordial 
non-Gaussianity induces a strong scale-dependence in the power and bi-spectra of biased objects 



at large length scales and changes the abun dance of very massive clusters (e.g. iDalal et al 



200S 



Nishimichi et al. 



20101 : 



Shandera et al. 



2011 



It is important to note that the large-scale clustering 
of galaxies and galaxy cluster can be used to constra in or possibly detect sc ale-dependent non- 



Gauss ianity which evades the CMB constraints (e.g. iLo Verde et al.l l2008l ). Ide Bernardis et al 



(j2010l ) derived constraints on /nl by combining the CMB anisotropies and galaxy clustering 



Croton et al. 


2007: 


Reid et al. 


2010a) 



galaxy bias, it would be ideal probe the dark matter distribution directly. 



Gravitational lensing is a powerful method to study dark matter distribution (e.g. lOguri et al 



2013). Future weak lensing surveys are aimed at measuring cosmic shear over a wide area of more 
than a thousand square degrees. Such observational programmes include Subaru Hyper Suprime- 
Cam (HSC) ^, the Dark Energy Survey (DES) ^, and the Large Synoptic Survey Telescope (LSST) 
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^. Space missions such as Euclid and WFIRST are also promising to conduct a very wide- field cos- 
mology survey. The large set of cosmic shear data will enable us to greatly improve the constraint s 



on cosmological parameters which include primordial non-Gaussianity (e.g. lOguri k Takadal[201ll ) 



There are many statistics proposed to characterize the large-scale matter distribution. Minkowski 
Functionals (MFs) are among useful statistics to extract the non-Gaussian information from two- 
dimensional or three-dimensional maps. For example, the full set of CM B MFs has alread y 
given comparable constr aints to those obtaine d using the CMB bispectrum (jHikage et al.ll2008l ). 
Matsubara fc JainI (120011) a nd ISato et al.l (|200ll ) studied ^rn-dependence of weak lensing MFs. More 
recently, Kratochvil et al. ( 20121 ) showed that the lensing MFs contain significant cosmological in- 
formation, beyond the power-spectrum. It is important and timely to study weak lensing MFs 
using fully nonlinear simulations of cosmic structure formation. 

In this paper, we forecast for future weak lensing surveys using MFs. In particular, we examine 
their ability to constrain the amplitude of the local-type primordial non-Gaussianity. We run a 
large set of A^-body simulations and then perform ray-tracing simulations of gravitational lensing. 
We measure MFs directly from the highly-resolved mock cosmic shear maps. We perform a full 
Fisher analysis to study the degeneracy among other cosmological parameters, especially the dark 
energy equation of state w = P/p. We propose to use future cosmology surveys to constrain, or 
possibly detect, primordial non-Gaussianities. 

The rest of the present paper is organized as follows. In Section [2l we describe the details 
of A^-body simulations and our ray-tracing simulations of gravitational lensing. In Section [3l we 
summarize the basics of MFs. In Section [H we show the results of a Fisher analysis using MFs. 
We clarify the degeneracy among the three parameters we consider. Concluding remarks and 
discussions are given in Section [5l 



METHODOLOGY 



2.1. iV-body Simulations 



We run a number of cosmological A'^-body simulations to generate w eak lensing convergence 
maps. We use the parallel Tree-Particle Mesh code Gadget 2 (ISpringelll2005l ) . Each simulation is run 
with 256^ dark matter particles. We run simulations of two different volumes, 240 and 480 /t~^Mpc 



on a s ide. We generate the initial conditi ons following a parallel code developed in iNishimichi et al 



([20091) and lValageas &: Nishiinichil (j201ll ). which employs the second-order Lagrangian perturbation 
theory (e.g. lCrocce et al.ll2006l ). The initial redshift is set to Zinit = 50, where we compute the linear 



^http://www.naoj.org/Projects/HSC/j_mdex.html 
^http: / /www. darkenergysurvey.org/ 
■^http: / /www. lsst.org/lsst / 



- 4 - 



matter transfer function using CAMB ([Lewis et alj |2000| ) . We then follow iNishimichi et al.l (|20ld ) 
to add non-Gaussian corrections to the initial conditions. For our fiducial cosmology, we adopt 
the following parameters: matter density ilm = 0.2726, dark energy density Q\ = 0.7274, Hubble 
parameter h = 0.704 and th e scalar spectral index = 0.963. These parameters are consistent with 
the WMAP 7- year results (iKomatsu et al.ll201ll ). For the primordial non-Gaussianity parameter, 
we adopt /nl = 0, ±100. To investigate the degeneracy of the cosmological parameters, we also 
run the same set of simulations for different w and erg, where as is the rms of the density field on 
8 /i~^Mpc. For these runs, we fix the amplitude of curvature fluctuations = 2.441 x 10~^ at 
the pivot scale k = 0.002Mpc~^. For w = —0.8, —1.0, and —1.2, the resulting erg is equal to 0.753, 
0.809, and 0.848, respectively. We summarize the simulation parameters in Tabled) 

After performing the simulations with w = —0.8 and w = —1.2, we found that our code 
had a bug in converting the physical time to the cosmic expansion parameter in the part where 
gravitational acceleration due to particle-particle interactions is calculated. This bug affected the 
results for w = —0.8 and w = —1.2. However, we have explicitly checked that the effect was 
very minor and that the statistics we use below were hardly affected. For example, the matter 
power spectra at z = after and before the bug was fixed differ less than 0.1 percent in the power 
amplitude at 0.04 < A: < 1 h/Mpc. 



2.2. Ray Tracing Simulations 



For ray-tracing simulations of gravitational lensing, we generate light-cone outputs using mul- 
tiple simulation boxes in the following manner. The small- and large- volume simulations are placed 
to cover a past light-cone of a hypothet ical observer with a ngula r extent 5° x 5°, from reds hift z = 
to z ^ 3.5, sim i larly to the methods in lwhite h Hul (|20od ) and iHamana h Mellierl (j200ll ). We fol- 
low lSato_et_alJ ([20091) in order to simulate gravitational lensing signals. Details of the configuration 
are found there. 

We set the initial ray directions on 2048^ grids. The corresponding angular grid size is 
5°/2048 ~ 0.15 arcmin. To avoid the same structure aligned along the line of sight, we shift 
randomly the A^-body simulation boxes. In addition, we use simulation outputs from independent 
realizations when generating the light-cone outputs. We generate 100 independent convergence 
maps from 20 A^-body simulations for each cosmological model. We fix the redshift of the source 

galaxies to 2:source = 

1.0,1.5. 

It is well-known that the intrinsic ellipticities of source galaxies induce noises to lensing shear 
maps. Assuming intrinsic ellipticities are uncorrelated, we compute the noise to convergence as 

('^noise(2^5 y)'^^noisc(^ i?/ )) — 7 ^xx'^yy'i (2) 

'^gal^pix 

where 5xx' is the Kronecker delta symbol, ngai is the number density of source g8ilB,xi6S , -^pix 
the solid angle of a pixel, and is the rms of the shear noise. Throughout this paper, we adopt 
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a-y = 0.4 and assume that the total number density of the source galaxies i s 30 galaxies /arcmin^ . 
These are typical values for a weak lensing survey using Subaru telescope (e.g. lMivazaki et al.ll2007l ). 
When we study a tomographic method (see Section [5]), we assume ngai = 15 galaxies/arcmin^ at 

-^source — 

1.0 and 1.5 0. To simulate a more realistic survey, we add the Gaussian noises following 
Eq.([2]) to our simulated maps. Then we perform the Gaussian filtering to the noisy lensing maps. 
We set the smoothing scale to 1 arcmin. This choice corresponds to th e optimal smoothing scale 
for the detection of massive halos using wGcik iGnsing with. -S^source — 

1.0 (|Hamana et al.ll2n04l V We 

discuss the effect of smoothing on the statistical analysis in Section O 



MINKOWSKI FUNCTIONALS 



3.1. Basics 



Minkowski Functionals are morphological statistics for some smoothed random field above a 
certain threshold. In general, for a given D-dimensional smoothed field, one can calculate D + 1 
MFs Vi. On S^, one can define 2+1 MFs Vb,Vi and V2. Vq, Vi and V2 describe the fraction of 
area above the threshold, the total boundary length of contours, and the integral of the geodesic 
curvature K along the contours. Mathematically, for a given threshold v, MFs are defined as 



1 
1 
1 



1 



-d£, 



1 

2^ 



KM, 



(3) 
(4) 
(5) 



where Qy and dQy represent the excursion set and the boundary of the excursion set for a smoothed 
field u{0). They are given by 



= {e\u{G) > u}, 

dQy = {e\u{e) = i^}. 



(6) 
(7) 



We follow iLim fc SimonI (j2012l ) to calculate the MFs from pixelated convergence maps. In this 
step, we convert a convergence field k into x = {k — (k)) / where (Tq is the standard deviation of a 
noisy convergence field on a 5° x 5° map. In binning the thresholds, we set Ax = 0.1 from x = —5 
to X = 5. We have checked that the binning is sufficient to reproduce the analytic MFs formula 



(jTomitalll986l ) for mock 1000 maps of Gaussian random fields. Figure [T] shows the measured and 
averaged MFs for our 100 convergence maps for the ACDM model and for ^source — 

1. We also 



With this simple split, the total noise per pixel can be kept constant. We can then study the significance of the 
high-z source galaxies for parameter constraints. 
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plot the analytic formula of MFs for Gaussian statistics to show the non-Gaussian features of the 
simulated convergence maps. 



3.2. Dependence on /] 



NL 



Let us first discuss how the primordial non-Gaussianity /nl affects the lensing convergence 
and the MFs. We define the ratio of MFs with respect to the fiducial /nl = model as follows: 



NLj 



V^ix■, /nl/0) 
V^ix■, /nL = 0)^ 



0,1,2. 



(8) 



Figure [2] shows Ri from our 100 convergence maps with ^source = 1-0. The effect of /nl appears 
large in the regime where the normalized convergence x > 3 for all ViS. This simply reflects the fact 
that /nl affects the nurn ber of very massive halos with mass ~ lO^^/i-^M© which yield x> 3 (see 
also iHamana et al.l |2004| ) . The abun dance of the massive h alos at z = 0.5 is larger by ~ 0(10%) 
for /nl = 100 compared to /nl = (INishimichi et al.ll2010l ). Then the fraction of area (Vq) with 
very high convergence increases for positive /nl, and the total length of contours (Vi) increases too. 
Interestingly, /nl also affects Vi and V2 at small x. Because the MFs, Vq, Vi, V2 are not independent 
statistics, their correlations need to be considered. We use all the MFs combined together in our 
statistical analysis below, in order to extract the full cosmological information and to derive an 
accurate constraint on /nl- 



4. RESULT 



4.1. Fisher analysis 

We perform a Fisher analysis to make forecasts for parameter constraints on /nL; w, and as 
for future weak lensing surveys. 

For a multivariate Gaussian likelihood, the Fisher matrix Fij can be written as 

F^j = ^Tr [AiAj + C-^Mij] , (9) 

where Ai = C^^dC /dS,i, Mij = 2 [d / d£,i) {d ^ / d^j) , C is the data covariance matrix, ^ is the 
assumed model, and ^ = (/nl, w-, cts)- For lensing MFs, /i corresponds to Vq, Vi, and V2 for 
different bins, n 



^We only consider the second term in Eq. ([9}. Since C is e xpected to scale i nverse-proportionally to the survey 
area, the second term will be dominant for a large area survey ijEifler et al.ll2009ll . 
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We estimate /i by averaging MFs over our 100 (noisy) convergence maps. To calculate Mjj, 
we approximate the first derivative of MFs by the cosmological parameter as follows 

We use the data set of the three MFs for ^source = 1-0 or/and 1.5. We use 10 bins in the range of 
X = [—3,3]. El We have checked this binning is sufficient to produce robust results in the following 
analysis. In this range of x, Eq. (jlOp gives smooth estimate for Mij. In total, we need 60 x 60 
MFs cov ariance matrix for the Fisher analysis. For this purpose, we use 1000 convergence maps 



made bv I Sato et al.l (j2009l ). These maps have the same design as our convergence maps, but are 



generated for slight ly different cosmological parameters (consistent with WMAP 3-years results 



(jSpergel et al.l 120071 )). We essentially assume that the dependence of the covariance matrix to 



cosmological parameters is unimportant. 

We also take into account the constraints from the CMB priors expected from the Planck 
satellite mission. When we compute the Fisher matrix for the CMB, we use the Markov-Chain 



Mont e-Carlo (MCMC) engine for exploring cosmological parameter space COSMOMC (jLewis &: Bridle 



20021 ). We consider the parameter constraints from the angular power spectra of temperature 
anisotropies, £'-mode polarization and their cross-correlation. For MCMC, in addition to ug and 
w, our independent variables include the matter density i^m^^) the baryon density Qbh'^, Hubble 
parameter h, reionization optical depth r, and the scalar spectral index Ug. To examine the pure 
power of lensing MFs to constrain /nL; we do not include any constraints on /nl from the CMB. 
Assuming that the constraints from the CMB and the lensing MFs are independent, we express the 
total Fisher matrix as 



MFs 



+ -f'cMB- (11) 



When we include the CMB priors by Eq. (jlip . we marginalize over the other cosmological param- 
eters except /nl; w and ag. 

Strictly speaking, one needs to consider a multivariate non-Gaussian likelihood because MFs 
are non-Gaussian estimator. In the present paper, we employ the Fisher analysis that assumes a 
local Gaussian likelihood in the parameter space. Note however that, because we use fully non- 
linear covariance matrices evaluated from 1000 ray-tracing simulations, our analysis appropriately 
includes non-Gaussian error contributions. Non-Gaussian error will include the contribution of 



four-point statistics at least (cf. iMunshi et al.ll2012bl ). It is illustrative to show the impact of non- 



Gaussian errors of MFs for parameter estimation. For comparison, we generate Gaussian covariance 
matrices by using 1000 Gaussian convergence maps. We have found that the resulting constraint 



^In principle, one could use regions with a; > 3, which are thought to be sensitive to /nl- However, such regions 
are extremely rare, and thus estimates for the first derivatives in Eq. (|10[) become uncertain even with our large 
number of convergence maps. 
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on cosmological parameters is degraded by a factor of a few percent compared to the case with 
Gaussian errors. 

4.2. Forecasts 

We show the forecast for the upcoming survey such as Subaru Hyper Suprime-Cam (HSC) 
and the Large Synoptic Survey Telescope (LSST). We consider two surveys with an area coverage 
of 1500 deg^ and 20000 deg^; the former is for HSC, and the latter is for LSST. We first derive 
constraints on the cosmological parameters for a 25 deg^ area survey, for which we have the full 
covariance matrix. Then we simply scale the covariance matrix by a factor of 25/1500 = 1/60 or 
25/20000 = 1/800 for the two surveys considered. 

Figure [3] shows the two-dimensional confidence contours for HSC (1500 deg^), in each case 
marginalized over other parameters. In each panel, the blue line shows the constraint from lensing 
MFs only and the red one is for the case of lensing MFs and the Planck priors. The ellipses shown 
in this figure correspond to the 68% confidence level from the Fisher analysis. Our fiducial model 
in this analysis is (/nl , w , as) = (0 , — 1 , 0.809). With the Planck priors, we can constrain /nl with 
a level of ~ 80 by using the source plane at ^source = 1-0, after marginalized over w and as- 

For the upcoming multiple-band imaging surveys such as HSC and LSST, it is possible to 
obtain photometric redshifts for the source galaxies. It is interesting to study how the parameter 
constraints can be improved by using source galaxies at higher redshifts. To this end, we perform 
the same analysis assuming the source galaxies are located at two redshifts, 2:source = 1-0 and 
1.5. We show the Fisher analysis result in the right panel of Figure El We see the constraint on 
/nl is improved by a factor of ~ 1.5 with the two redshift bins. We summarize the marginalized 
constraints on /nl in Table [2j 

5. CONCLUSION AND DISCUSSION 

We have studied the ability of weak lensing MFs to constrain the local type primordial non- 
Gaussianity. We have performed 20 A^-body simulations and 100 independent ray-tracing sim- 
ulations for each of the seven cosmological models that differ in /nl, w and as- We have then 
performed a Fisher analysis using the large set of mock lensing maps, to obtain confidence limits 
for cosmological parameters. 

The MFs are sensitive probes of /nl, especially at high convergence values, because such high 
convergence regions are associated with massive halos with > 10^^h~^ Mq, of which the abundance 
is sensitive to /nl- We also find that the three MFs, Vq, Vi and V2, are affected differently by /nl- 
This means that combining the three MFs gives tighter constraints on /nl- 

From a Fisher analysis, we have obtained the following results. For source galaxies at ^source = 



- 9 - 



1.0, the primordial non-Gaussianity is constrained with a level of /nl ~ 80 for HSC survey with a 
1500 deg^ survey area. The constraints can be improved by selecting source galaxies in multiple 
redshift bins. We find that the constraint on /nl is improved to ~ 50, i.e., by a factor of ~ 1.5 if we 
use source galaxies at ^source = 1-0 and 1.5. This is largely because the degeneracy between as and 
w is broken due to information contained in the matter distribution at different redshifts. We have 
also tested how the number density of the source galaxies affects our analysis. We have re-analyzed 
the case of rigai = 30 galaxies/arcmin^ for a fixed source plane at ^^source = 1-0. The constraint 
on /nl is improved only by a factor of ~ 5 % in this case. We argue that the "tomographic" 
information using the multiple source planes is useful to derive accurate constraints on /nl, even 
though one then needs to use a smaller number of source galaxies at each of the source planes. 

Ultimately, for a LSST-like survey with a 20000 deg^ area, we can obtain the constraint of 
/nl ~ 25 with ^source = 1-0, and /nl ~ 15 with Zsource = 1-0 and 1.5. In principle, these constraints 
will be further improved by including the high sigma bins of MFs because the higher convergence 
region is more sensitive to /nl- We will continue our study along this line using a larger set of 
simulations of larger volumes. 

Finally, we discuss possible technical improvements in using the MFs of weak lensing maps. 
First, we have checked how the smoothing scale on the convergence maps affects our analysis. We 
have performed a Fisher analysis using the maps with smoothing of 0.5 arcmin, 1 arcmin, 2 arcmin 
and 5 arcmin. Smoothing with ~ 1 arcmin has turned out to be optimal for the targeted future 
surveys, yielding the best constraints. This is explained qualitatively by the fact that only the large- 
scale, linear structure is probed with large smoothing scales, whereas for smaller smoothing scales, 
the intrinsic shape noise becomes large. We note that the angular-scale dependence of MFs itself 
can provide more information. For example, the angular-scale dependence can be used to separa te 
primordial non-Gaussi anities from gravity-induced non-Gaussianities (e.g. Munshi et al. 2012al lbl: 



Kratochvil et al.l |2012| ) . One could further improve the cosmological constraints by combining 
MFs with various smoothing scales and their evolutions. Evaluating MFs is complicated when an 
observed map include masked regions. Also, instrumental and atmospheric systematics can easily 
compromise the measurement of lensing MFs. These issues certainly warrants further extensive 
studies. The upcoming wide-field surveys will provide highly-resolved lensing maps. Our study in 
the present paper may be useful to properly analyze the data and extract cosmological information 
from them. 
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Table 1: Parameters for our A^-body simulations. For each model, we run 20 A^-body reahzations 
and generate 100 weak lensing convergence maps. 
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Fig. 1. — The Minkowski Functionals from our simulated convergence maps. The solid lines show 
the average MFs over 100 maps with the fiducial co smological pa rameters. The dashed hues show 
the analytic formula of MFs for the Gaussian field (lTomitalll986l ). We use 100 convergence maps 
located on 2;soiirrp = 1- 





^source — I'O 


^source — 1.0, 1.5 


MFs only (1500 deg^) 


96.8 


57.9 


MFs + Planck (1500 deg^) 


78.5 


52.1 


MFs only (20000 deg^) 


26.5 


15.8 


MFs + Planck (20000 deg^) 


25.8 


15.7 



Table 2: The 1-cr constraint on /nl when marginalized over w and erg. We consider two surveys 
with a survey area of 1500 deg^ (HSC) and 20000 deg^ (LSST). The analysis includes the intrinsic 
noise from source galaxies with the number density of ngai = 15 galaxies/arcmin^ at each source 
redshift. 



- 13 - 




Fig. 2. — The MFs for our simulated maps with /nl- The vertical axis shows the ratio of lensing 
MFs 1^(/nl = ilOO)/T^(/NL = 0). The horizontal axis shows the normalized convergence field 
(k — (k))/cjo. The red (blue) points with error bar represent the result from 100 noisy convergence 
maps with /nl = +100 (-100). The source plane locates on ^source = 1-0 and 1 arcmin Gaussian 
smoothing is adopted. 
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Fig. 3. — We plot 1-cr confidence level by weak lensing MFs for Subaru HSC survey (1500 deg^). 
The left panel shows the constraints on /nl, w and erg for the case with ^source = 1-0. The right 
panel shows the constraints obtained by a tomographic survey with ^source = 1-0 and 1.5. The blue 
line shows the constraints from lensing MFs only and the red one represents those by MFs and the 
Planck priors. 



